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Abstract 

I present results for top quark production in hadronic collisions at LHC and Tevatron 
energies. The soft-gluon corrections to the differential cross section are resummed at next- 
to-next-to-leading-logarithm (NNLL) accuracy via the two-loop soft anomalous dimension 
matrices. Approximate next-to-next-to-leading-order (NNLO) differential and total cross 
sections are calculated. Detailed theoretical predictions are shown for the ti cross section 



Q H i and the top quark px distribution at the Tevatron and the LHC. 
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J . 1 Introduction 

The top quark occupies a unique position in the list of elementary particles as the most massive 
particle discovered to date. Its high mass suggests an important role for the top quark in the 
physics of electroweak symmetry breaking. After a long period of searches, the discovery of the 
top quark via top-antitop production in proton-antiproton collisions (pp — > ti) was announced 
in 1995 by the CDF and DO collaborations at the Fermilab Tevatron collider pp. The ti cross 



On 

section has been measured with increasing precision at Run II of the Tevatron [21 [3] and there 
has also been data for the transverse momentum, px, distribution of the top quark [4]. More 
recently single top quark production was observed by DO [5] and CDF [6] . Measurements of the 
top quark mass have also been increasingly more precise [7j. The LHC is expected to observe a 
very large number of top quark events and to bring top quark physics to a new energy frontier. 
For reviews of top quark physics at the Tevatron and the LHC see Ref. [8] (experiment) and 
Ref. [9] (theory). 

The experimental measurements of the top quark cross section and px distribution at the 
Tevatron are currently in good agreement with theoretical predictions. However, as the experi- 
mental errors continue to get smaller with time, precise theoretical calculations with smaller un- 
certainties are required. Next-to-leading order (NLO) calculations of the QCD corrections have 
been available for over two decades [T0~1 [TT] (electroweak corrections, which are much smaller 
numerically, have also been calculated more recently [12]) but the associated uncertainty is 
much bigger than current experimental errors. The inclusion of higher-order soft-gluon correc- 
tions enhances the cross section and px distribution and significantly reduces the theoretical 
error [HUH]. 

Until recently, the state of the art in theoretical predictions was approximate next-to-next- 
to-leading order (NNLO) calculations based on next-to-leading-logarithm (NLL) resummation 
of soft-gluon corrections for the differential cross section, supplemented with further subleading 



1 



terms [T3| [T3] . These soft-gluon corrections are dominant not only near partonic threshold but 
also away from it. The accuracy at NLL was achieved by the calculation of the one- loop soft 
anomalous dimension matrices for the partonic channels in top quark production in Ref. [15]. 

To achieve next-to-next-to-leading-logarithm (NNLL) accuracy in the resummation one 
needs to calculate the soft anomalous dimensions at two loops. This is a much more difficult 
undertaking. For massless quark scattering, the two-loop soft anomalous dimension matrix was 
first calculated in Ref. [16]. Further work on soft and collinear singularities of dimensionally- 
regularized scattering amplitudes in massless gauge theories followed in Refs. [T7J HHJ [191 I2Q1 [21] . 
More recently, a lot of work on massive two-loop soft anomalous dimensions has appeared in 
Refs. [221 EHl EH EHl EEl [271 EHl EH1 EU] EU [32]. The presence of a mass for the top quark 
considerably complicates the calculation relative to the massless case. 

Soft-gluon resummation is a consequence of factorization. The partonic cross section can be 
factorized into functions associated with the hard scattering, collinear and soft-gluon emission 
from the external partons, and noncollinear soft-gluon emission that depends on the color 
structure of the process [To]. The renormalization group evolution of these functions results 
in expressions for the resummed cross section. The resummation formalism followed here has 
already been presented and reviewed in numerous papers over a period of more than a decade 
(see Refs. [131 US [23 [291 ESI [Ml [35] and references therein) so we will not repeat the derivation 
of resummation and we will not repeat explicit expressions in this paper except for new two- 
loop results in Section 2. Resummation is performed in Mellin moment space: we define a 
kinematical variable that measures distance from partonic threshold, and then iV is the 
moment variable conjugate to S4. For ti production the resummed partonic cross section in 
moment space is given by 



The first exponent in the above expression resums soft and collinear corrections from the in- 
coming partons a and b (quark-antiquark or gluon-gluon) while the second exponent controls 
the factorization scale, fip, dependence of the cross section. H ab is the hard-scattering function 
while S a b is the soft function describing noncollinear soft gluon emission. The renormalization 
group evolution of the soft function is controlled by the soft anomalous dimension, T Sa b [15J. 
It is important to note that H a b, S a b, and Ts a b are matrices in the space of color structures 
of the process [T5l [331 EI]- I n the next section we will present explicit expressions for the new 
two-loop results for the soft anomalous dimension matrices ^sqg, for the qq — > ti channel, and 
^Sgg, for the gg — > ti channel. It is these new ingredients that allow us to complete the NNLL 
resummation in our formalism (for other approaches see Refs. [251 ESI E3 EH [32] and the 
discussion in Section 5). 

The resummed cross section, Eq. (jl.ip can be expanded at fixed order in a s to NLO, NNLO, 
etc., and inverted back to momentum space, see e.g. Refs. [131 EH EH [35] . At each order in a s , 
one encounters plus-distribution terms of the form [ln fc (s4/m 2 )/s 4 ] + , where m is the top quark 
mass and, for the n-th order corrections, the power of the logarithm, k, can range from the 
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leading value of 2n — 1 down to the lowest value of 0. Thus at NLO k = 1 or 0, while at NNLO 
k can take the values 3,2,1,0. From NLL resummation one can determine the coefficients of 
both, k — 1,0, powers of the logarithms at NLO, but only the powers k — 3, 2, 1 at NNLO 
(determining the NNLO k = 1 term requires matching with NLO). Partial results for the k = 
term at NNLO were provided in Ref. [13] and also used in [13]. From NNLL resummation one 
can in addition fully determine the k = term at NNLO. 

In the following section we present the soft anomalous dimension matrices for the qq — > ti 
and gg — > ti channels at one and two loops. In Section 3 we use the NNLL resummation to 
obtain approximate NNLO results for the total ti cross section and the top quark px distribution 
in proton-antiproton collisions at the Tevatron. In Section 4 corresponding results are given for 
proton-proton collisions at LHC energies. A comparison with other approaches and conclusions 
are given in Section 5. 



2 Soft anomalous dimension matrices for tt production 



We begin with the result for the soft (cusp) anomalous dimension T$ [22\ for e + e~ — > ti, which is 
an integral part of the calculation for the soft anomalous dimension matrices Ts q q and ^sgg f° r ti 
hadroproduction. The calculations of soft anomalous dimensions involve diagrams with eikonal 
lines representing the top quarks. The eikonal diagrams are calculated in Feynman gauge in 
momentum space, and we use dimensional regular izat ion with d = 4 — e dimensions to isolate 
the ultraviolet (UV) poles of the diagrams. The soft anomalous dimension is then determined 
from the coefficients of the UV poles [22] • Writing T s = (a s /iT)T^ + (a s /7r) 2 Tg H — •, we have 
the one-loop expression 
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The two-loop soft (cusp) anomalous dimension, determined from the UV poles of two-loop 
eikonal diagrams, is [221 121] 
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where K = (7^(67/18 — (2) — 5ri//9, with Ca = N c and n/ = 5 the number of light-quark 
flavors. We have written the two-loop result in Eq. (12. 3p in the form of a term which is a 
multiple of the one- loop soft anomalous dimension ri 1 ^ plus additional terms, denoted as Ma: 
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This result, first obtained in [22], is written in terms of logarithms, dilogarithms, and triloga- 
rithms, and it provides a more explicit analytical expression than earlier work [36J. Note that 
as $ 1, M p (1 - Cs)/2. 

We can now proceed with the results for the two-loop soft anomalous dimension matrices 
for the partonic processes qq — > tt and gg — >■ tt. The calculation involves the two-loop soft 
(cusp) anomalous dimension for all pairs of external lines in the process (cf. [221 ESI 29J) as 
well as graphs with gluons connecting three external lines (cf. (23J, [27J [31]). We begin with top 
quark production through light quark annihilation, 



q(Pa) + q(pb) -> t(pi) + t(p 2 ) ■ 
We define the kinematical invariants 

S = (p a +Pb) 2 , h = (Pb-Pl) 2 -m 2 , Ul = (p a 
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and S4 = s + ti + Ui, where S4 measures distance from partonic threshold. The calculations are 
performed in a color tensor basis consisting of singlet and octet exchange in the s channel, 
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Here the color indices for the incoming (light) quark and antiquark are a and b, respectively, 
and for the outgoing top quark and antiquark 1 and 2, respectively, and Tp are the generators 
of SU(3) in the fundamental representation. 

The matrix for qq — > tt in this c\, c 2 color basis is 
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The result in Eq. ( 12. 9p is somewhat different from the original in Ref. [15] because the 
original calculation used the axial gauge while Eq. ( 12. 9 p is in Feynman gauge. Of course this 
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does not affect the complete resummed expression because other terms in the resummed cross 
section compensate by also taking different forms in the two gauges. We note that the "11" 
element of the matrix is simply the cusp anomalous dimension, Tg. 
At two loops: 
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comes from graphs with gluons connecting three external lines, whose contribution were first 
calculated explicitly in [27]. Note that Np is just a subset of the terms of Mp, Eq. ( 12. 4p . so all 
analytical structures already appear in Mp, and that as /3 — > 1, iVg — > 0. The two-loop matrix, 
Eq. 1 12. 10p is not proportional to the one-loop matrix, Eq. 1 12. 91) . This fact was first discussed 
in Ref. |22j and it is to be contrasted with the simple proportionality relation for the massless 
case that was found in Ref. [T6] . 
We continue with the gg channel: 



g(p a ) + g{Pb) -> *(pi) + <(P2) • 
We choose the following basis for the color factors: 
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where d abc and / abc are the totally symmetric and antisymmetric SU(3) invariant tensors, 
respectively. We define s, t X) and u\ for this channel as in Eq. 112.61) . 
The matrix for gg — > ti in this basis is 
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The expression in Eq. ( 12 .151) is again somewhat different from the original in Ref. [15] because 
Eq. (I2.15P is derived in Feynman gauge. 
At two loops: 
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As was the case for the channel, we note that for the gg channel the two-loop matrix, Eq. 
(12.161) is not proportional to the one-loop matrix, Eq. (I2.15p . 

The expressions in Eqs. (12.101) and (I2.16P are different from the corresponding ones in [27] 
due to different definitions and formalism. 

With the two-loop soft anomalous dimension matrices at hand we achieve NNLL accuracy 
in the resummed cross section, Eq. (11.11) . Expanding the resummed cross section to NNLO we 
then calculate approximate NNLO cross sections and transverse momentum distributions for 
top quarks at the Tevatron and the LHC. 



3 Top cross section and £>t distribution at the Tevatron 

We now provide a detailed phenomenological study of top quark production at the Tevatron 
collider, including the total ti cross section and the top quark pt distribution. We present NLO 
and approximate NNLO calculations for these quantities. The NNLO approximate results are 
computed by adding the NNLO soft-gluon corrections (derived from NNLL resummation) to 
the exact NLO quantities. The total and differential cross sections depend on the factorization 
scale, fMp, and the renormalization scale, fiR. These two scales are often set equal to each other 
and denoted simply as /x, but they are in principle independent. 
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Figure 1: The NLO and approximate NNLO cross section for tt production at the Tevatron 
with VS = 1.96 TeV and MSTW2008 NNLO pdf. 



3.1 tt cross section at the Tevatron 

In Fig. [TJ we plot the cross section for top-antitop production in proton-antiproton collisions 
at the Tevatron over a top quark mass range 165 < m < 180 GeV at a factorization and 
renormalization scale fi = m. The exact NLO and the approximate NNLO cross sections are 
shown. The enhancement from the NNLO soft-gluon corrections is 7.8%. Here we have used 
the MSTW2008 NNLO parton distibution functions (pdf) [37]. We will use these pdf for our 
calculations throughout this paper except where noted otherwise. 

Table 1 lists the values for the NNLO approximate cross section at the Tevatron for top 
quark masses between 170 GeV and 175 GeV. There are two kinds of theoretical uncertainties 
associated with the calculation: dependence on the factorization/renormalization scale, and 
uncertainties from the parton densities. 

The scale dependence of the cross section for m = 173 GeV is plotted in Fig. [2] over a 
range of two orders of magnitude, 0.2 < [i/m < 10. It is clear that at leading order (LO) 
the cross section is strongly dependent on the choice of scale, varying by a factor of 5.45 
between maximum and minimum values in the range shown. The NLO corrections significantly 
stabilize the LO variation: the NLO cross section varies by a factor of 1.61. The NNLO soft- 
gluon corrections further reduce the scale dependence: the NNLO approximate cross section 
varies by a factor of only 1.18. The improvement provided by the NNLO corrections is even 
more impressive if one considers only the variation 0.5 < \xjm < 2 as traditionally used to 
estimate errors. For this range the LO cross section varies by a factor of 1.85, the NLO cross 
section by 1.16, while the NNLO approximate cross section by a factor of only 1.034. 
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NNLO approx tt cross section (pb) 


m (GeV) 


Tevatron 


LHC 7 TeV 


LHC 14 TeV 


170 


7.78 


179 


998 


171 


7.54 


173 


972 


172 


7.31 


168 


946 


173 


7.08 


163 


920 


174 


6.87 


158 


896 


175 


6.66 


154 


873 



Table 1: The NNLO approximate tt production cross section in pb in pp collisions at the 
Tevatron with =1.96 TeV and in pp collisions at the LHC with \/S = 7 TeV and 14 TeV. 
We set y. = m and use the MSTW2008 NNLO pdf [37]. 



For a top quark mass of 173 GeV, the NLO cross section is 6.57l ; 66 _ ; 25 pb and the NNLO 
approximate cross section is 

a NNLOapprox (m = m ^ j g g ^ = 7 .Qg^f^fr pb . (3.1) 

Here the first uncertainty is from scale variation over 0.5 < fi/m < 2 and the second is 
from the MSTW2008 NNLO pdf errors at 90% C.L. (to be conservative, we do not use the 
smaller 68% C.L. pdf errors). At NLO the scale uncertainty is bigger than that from the 
pdf, but at NNLO the scale uncertainty is much smaller than the pdf one. In fact the scale 
uncertainty at NNLO is about four times smaller than that at NLO, again highlighting the 
dramatic reduction of scale dependence provided by the higher-order corrections. Adding the 
scale and pdf errors in quadrature, the NNLO approximate result is 7.08 ±0.36 pb, i.e. we have 
a ± 5.1% total uncertainty, which is to be contrasted with a much larger (+6.6% -10.7%) total 
error (in quadrature) at NLO. 

One can also study the dependence of the cross section separately on the factorization scale 
and the renormalization scale. This can be important because in some cases setting \xp equal 
to \xr may give a smaller uncertainty than from varying the scales independently. In Fig. [3] we 
plot the scale dependence of the cross section for m = 173 GeV in three different ways at NLO 
(top plot) and approximate NNLO (bottom plot). The first way is to set /i = fip = and vary 
this common scale, exactly as we did in Fig. [2J The second way is to vary the factorization scale 
fip while keeping the renormalization scale fixed at /j,r = m. The third way is to vary fiR while 
keeping \ip — m. It is clear from the top plot that varying fip and [ir independenty over the 
range m/2 and 2m does not give a wider range of cross section values than varying the common 
scale /i — i~Lf — I-Lr- In fact as can be seen from the figure this holds true for a very wide range 
of scale variation. We also note that setting fip = m/2 and fifi = 2m or setting fiR = m/2 and 
fip = 2m still gives a smaller variation than varying the common scale fi = fip = fiR between 
m/2 and 2m. Therefore the NLO theoretical uncertainty that we provided above from scale 
variation is not increased by separately varying /j,p and //r. For the approximate NNLO cross 
section in the bottom plot of Fig. [3] we see that the variation with /j, = \ip and \ir = m affects 
the upper uncertainty (which was stated before as +0.00) and this new upper uncertainty is 
+0.20. However the lower uncertainty (—0.24) is unaffected. So the result for the approximate 
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pp -> tt at Tevatron S =1.96 TeV m=173 GeV 
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Figure 2: The scale dependence of the tt cross section at the Tevatron with v^S* = 1.96 TeV 
and m = 173 GeV. 
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Figure 3: The fip and fiR dependence of the tt cross section at the Tevatron with vS* = 1.96 
TeV and m = 173 GeV. The top plot is at NLO and the bottom is at approximate NNLO 
accuracy. 
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pp -> tt at Tevatron S =1.96 TeV m=173 GeV 
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Figure 4: The top quark p T distribution at the Tevatron with y/S = 1.96 TeV and m = 173 
GeV. 



NNLO cross section for m = 173 GeV with the scale uncertainy from independent up and [1r 
variation can be written as 7.08i ]24io!27 pb- Finally, we note that not only is the scale variation 
with p, = fip = fiR greatly reduced in going from NLO to approximate NNLO but so is the 
separate /xp variation and the separate /Jr variation. 

The MSTW2008 parton densities are the only ones available at NNLO and so we use them 
for our best predictions. It is interesting nevertheless to see if the results change significantly 
using the new CT10 pdf [38J, which are at NLO, and the pdf errors associated with them. 
Using CT10 pdf we find a NLO cross section for m = 173 GeV of 6.81±g;fl±g;|§ pb, and an 
approximate NNLO cross section of 7.38^o;25^o;32 pb, where the first uncertainty is from scale 
variation (with /ip and [Ir independently varied) and the second is from the pdf errors. We thus 
find both a larger cross section and a larger uncertainty with CT10 pdf than with MSTW2008 
NNLO pdf. 

3.2 Top quark distribution at the Tevatron 

The top quark transverse momentum distribution at the Tevatron with m = 173 GeV is plotted 
in Figs. HI and using the MSTW2008 NNLO pdf. Fig. H shows the differential distribution 
da/dpT over a range < pt < 300 GeV. Both NLO and NNLO approximate results are shown 
for three different scale choices, /i = m/2, m, and 2m. The integrated pt distribution gives the 
same result for the total cross section as found in the previous subsection, which provides a 
good consistency check of the calculation. The scale variation of the pt distribution at NNLO 
is again significantly smaller than at NLO. The NNLO soft-gluon corrections enhance the NLO 
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pp -> tt at Tevatron S =1.96 TeV m=173 GeV 
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Figure 5: The top quark px distribution at the Tevatron with y/S = 1.96 TeV, m = 173 GeV, 
and \x = m or fi = in a logarithmic plot. 



result but the shape is similar. 

Figure [5] presents the top quark distribution in a logarithmic plot that makes it easier 
to see da/dpT at high p T values. Results are now shown up to a px of 500 GeV. In Fig. H]the 
central value for the scale was taken to be fi = m as for the total cross section, and the scale 
variation was around that central value. Another possible scale choice for the top quark p? 
distribution is the transverse mass mr, defined by mr = (pt + m 2 ) 1 / 2 . In Fig. [5] we show our 
NLO and approximate NNLO results for both fi = m and fi = tjit- We find that the choice of 
scale, m versus m^, makes very little difference even for high p? of 500 GeV - the curves are 
practically indistinguishable. 

Joint threshold and recoil resummation for the pt distribution (at NLL accuracy only) has 
been studied in [39J. The effect of recoil is entirely negligible except at extremely high pt (~800 
GeV and above) so we do not consider it further. 

4 Top cross section and distribution at the LHC 

We continue with a detailed phenomenological study of top quark production in proton-proton 
collisions at the LHC. We present results for the current LHC energy of 7 TeV and the future 
(design) energy of 14 TeV, and also a few results at 10 TeV. 
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pp -> tt at LHC S =7 TeV |a=m 



220 



200 - 



NNLO approx 




120 - 




170 



175 



180 



m (GeV) 



Figure 6: The cross section for tt production at the LHC with = 7 TeV and MSTW2008 
NNLO pdf. 



4.1 tt cross section at the LHC 

In Fig. [6] we plot the NLO and approximate NNLO cross section for top-antitop production at 
the LHC at 7 TeV energy over a top quark mass range 165 < m < 180 GeV at a factorization 
and renormalization scale // = m using the MSTW2008 NNLO pdf. The enhancement from 
the NNLO soft-gluon corrections is 7.6%. Table 1 lists the values for the NNLO approximate 
cross section at the LHC at an enery of 7 TeV for top quark masses between 170 GeV and 175 



The scale dependence of the cross section for m = 173 GeV is plotted in Fig. [7] over 
a range of two orders of magnitude, 0.2 < ///to < 10. Again, at LO the cross section is 
strongly dependent on the choice of scale, varying by a factor of 4.64 between maximum and 
minimum values in the range shown. The NLO corrections stabilize the LO variation: the 
NLO cross section varies by a factor of 1.85. The NNLO soft-gluon corrections further reduce 
the scale dependence: the NNLO approximate cross section varies by a factor of 1.43. The 
improvement from the NNLO corrections is again more impressive if one considers only the 
traditional variation 0.5 < fi/m < 2. For this range the LO cross section varies by a factor of 
1.75, the NLO cross section by 1.27, while the NNLO approximate cross section by a factor of 
only 1.08. 

For a top quark mass of 173 GeV, the NLO cross section is 152±i|±| pb and the NNLO 
approximate cross section is 



GeV. 



.NNLOapprox 
tt 



(to = 173 GeV, 7 TeV) = 163±^±| pb 



(4.1) 
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pp -> tt at LHC S =7 TeV m=173 GeV 
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Figure 7: The scale dependence of the tt cross section at the LHC with \fS = 7 TeV and 
m = 173 GeV. 



where the first uncertainty is from scale variation over 0.5 < fi/m < 2 and the second is from 
the MSTW2008 NNLO pdf errors at 90% C.L. At NLO the scale uncertainty is about twice as 
big as the pdf one, but at NNLO it is significantly smaller. The scale uncertainty at NNLO is 
about three times smaller than that at NLO. Adding the scale and pdf errors in quadrature the 
NNLO approximate result is 163j^o pb, i.e. we have a +7.0% -6.3% total uncertainty, which is 
to be contrasted with a much larger (+11.8% -13.8%) total error at NLO. 

We also study the dependence of the cross section separately on the factorization scale and 
the renormalization scale. In Fig. [8] we plot the scale dependence of the cross section in three 
different ways at NLO (top plot) and approximate NNLO (bottom plot). The first way is to 
set fi = \ip = //r and vary this common scale, as we did in Fig. [71 The second way is to 
vary fip while keeping \ir = m, and the third way is to vary fiR while keeping fip = m. From 
the top plot we see that varying fi F and fiR independenty over the range m/2 and 2m does 
not give a wider range of cross section values than varying the common scale fi = \ip = [J,r, 
and this actually holds true for nearly the entire wide range of scale variation shown in the 
plot. We also find that setting y,p = m/2 and Hr = 2m or setting fiR = m/2 and fip = 2m 
still gives a smaller variation than varying the common scale fi = fip = between m/2 and 
2m. Therefore the NLO theoretical uncertainty from scale variation provided previously is not 
increased by separately varying fip and fiR. For the approximate NNLO cross section in the 
bottom plot of Fig. [8] we also see that the independent variation of hf and hr does not affect 
the uncertainty that we wrote previously. Finally, we note that the separate hf variation and 
[ir variation are reduced when going from NLO to approximate NNLO. 

Again we can check if the results change significantly using the CT10 pdf, which are at 
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Figure 8: The [ip and hr dependence of the tt cross section at the LHC with \/lS = 7 TeV and 
m = 173 GeV. The top plot is at NLO and the bottom is at approximate NNLO accuracy. 



NLO. Using CT10 pdf we find a NLO cross section at m = 173 GeV of 15013,0^10 pb, an d an 
approximate NNLO cross section of 162^1]^ pb, so the results are very similar to those with 
MSTW2008 NNLO pdf. 

For reference, the cross section at a possible future LHC energy of 10 TeV is plotted in 
Fig. M using the MSTW2008 NNLO pdf. For a top quark mass of 173 GeV, we find a NLO 
cross section of 385^45l}g pb, while at NNLO 

a NNLOapprox (m = ^3 ^ 1Q ^ = 415 +17+18 pb ^ 

The cross section for the design LHC energy of 14 TeV is plotted in Fig. [10] using the 
MSTW2008 NNLO pdf. The enhancement from the NNLO soft-gluon corrections is 8.0%. 
Table 1 lists the values for the NNLO approximate cross section at 14 TeV LHC energy for top 
quark masses between 170 GeV and 175 GeV. The NLO cross section for a top quark mass of 
173 GeV is 852^1| pb and the approximate NNLO cross section is 

a NNLOapprox (m = ^3 ^ ^ ^ = 920+^+33 pb . (4.3) 

Again we observe a significant decrease in scale dependence at NNLO relative to NLO, and 
also note that separate variation of /ip and hr does not increase the uncertainty. The pdf 
uncertainties at this high energy are much smaller than the scale variation at NLO, and some- 
what relatively smaller at NNLO. Adding the scale and pdf errors in quadrature the NNLO 
approximate result is 920^52 pb, i.e. we have a +6.5% -5.7% total uncertainty, which is to be 
contrasted with a much larger (+11.2% -11.6%) total error at NLO. 
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Figure 9: The cross section for tt production at the LHC with = 10 TeV and MSTW200£ 
NNLO pdf. 
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Figure 10: The cross section for tt production at the LHC with y/S = 14 TeV and MSTW200£ 
NNLO pdf. 
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Figure 11: The top quark distribution at the LHC with v^S* = 7 TeV and m = 173 GeV. 



4.2 Top quark pr distribution at the LHC 

The transverse momentum distribution of the top quark with m = 173 GeV at the LHC at 
7 TeV energy is plotted in Figs. M and QJ using the MSTW2008 NNLO pdf. Fig. [TT] shows 
NLO and approximate NNLO results for the differential distribution da/dp? over a range 
< Pt < 350 GeV for three different scale choices, fi — m/2, m, and 2m. The scale variation 
of the pt distribution at NNLO is much smaller than that at NLO. 

Figure [12] presents the results for da /dp? in a logarithmic plot for high p? values up to 
1000 GeV, using both fi = m and fi = rriT, where again is the transverse mass. At very 
high pj- the NNLO approximate corrections become increasingly more significant and begin to 
change the shape of the distribution relative to NLO. This is not unexpected since the soft- 
gluon corrections are large near partonic threshold, which is dominant at high p?. The change 
of shape is more pronounced with the choice fi = m than it is with fi = ttlt- 

The pt distribution of the top quark with m = 173 GeV at the LHC at 14 TeV energy is 
plotted in Figs. [13] and [HI Fig. [13] shows NLO and approximate NNLO results over a range 
< Pt < 400 GeV for three different scale choices, fi = m/2, m, and 2m. Again, the scale 
variation of the pt distribution at NNLO is much smaller than that at NLO. 

Figure [TH presents the results for da/dpT in a logarithmic plot for high pt values up to 
1500 GeV, using fi = m and fi = mT- The NNLO soft-gluon corrections provide a significant 
enhancement and change the shape of the NLO distribution at very high pt- Again, the change 
in shape is larger with the choice fi = m than it is with fi = mT- 
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Figure 12: The top quark p T distribution at the LHC with y/S = 7 TeV, m = 173 GeV, and 
/i = m or |j = tut in a logarithmic plot. 
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Figure 13: The top quark px distribution at the LHC with v^S* = 14 TeV and m = 173 GeV. 
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Figure 14: The top quark p T distribution at the LHC with y/S = 14 TeV, m = 173 GeV, and 
fi = m or /i = rriT in a logarithmic plot. 

5 Comparison with other approaches and Conclusions 

In this paper we have resummed the soft-gluon logarithms in top quark production to NNLL 
accuracy. This work directly extends the earlier NLL resummation in Ref. [13] and the further 
work in [321 EH] and later in [T3"| IT4"] . To achieve NNLL (NLL) accuracy, it is necessary to 
derive the soft anomalous dimension matrix at two (one) loops for each partonic process. At 
NLL and NNLL the color structure of the hard scattering enters the resummation in a non- 
trivial way. The soft anomalous dimension matrices are explicitly dependent on the kinematical 
variables s, ti, Ui, and the resummation involves these quantities and logarithms of s^, where 
s 4 = s + ti + Ui measures distance from partonic threshold. Thus this is a fully differential cal- 
culation and the formalism in this paper has been used to calculate not only total cross sections 
but also differential cross sections, such as transverse momentum distributions. Approximate 
NNLO differential cross sections are extracted from the resummation (higher-order contribu- 
tions beyond NNLO are small, see e.g. Ref. [35J). The NNLO expansion avoids the need for 
prescriptions to deal with Landau-pole divergences in the resummation, and we prefer to take 
this approach since the numerical discrepancies between different prescriptions are larger than 
the corrections beyond NNLO (see e.g. discussion in [33| [14"]). 

There also exist formalisms of resummation and finite-order expansions for the total cross 
section only [ID] that are calculationally simpler, and only involve the variable (3 = yjl — 4m 2 /s. 
Logarithms of have been resummed at NLL in [10J HI] and at (partial) NNLL in [12] (Ref. 
[42J made an incorrect assumption about the two-loop terms which, as later understood, is not 
valid). This approach does not, however, involve the exact differential kinematics and hence 
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numerical deviations from the exact kinematics-sensitive result may appear. Furthermore this 
approach is inapplicable to pr or other differential distributions, so it is limited in scope. For 
further discussion of the differences see also Ref. [H] and [30]. More recently, complete NNLL 
results in logarithms of /3 have appeared in Ref. (25] using soft-collinear and non-relativistic 
effective theory, and in Ref. [2B] using resummation in moment space. Threshold expansions 
to NNLO for the total cross section from this ln/3 resummation have recently appeared in 
|l3l HI]. Again, all these results are for total cross sections only, based on expansions in (3. It 
is important to note that the terminology "NLL" and "NNLL" means different things in the 
approaches of Refs. [251 1261 HO, HD H21 H3J SI] than it does in the differential-level formalism of 
Refs. [131 HH [El E31 El] and of this paper because different types of logarithms are resummed. 

Another differential-level formalism that has recently appeared [30| |4"5] is based on soft- 
collinear effective theory and heavy-quark effective theory. While the resummation in [T3| [HI 
IT5j 1331 El] and this paper is done in moment space, in [30] it is performed in momentum space. 
The total cross section and invariant mass distribution at NNLL have been presented in [30J. 
The total cross section results in [30J are quite different from those in this paper. One major 
reason for the difference is the different choice of kinematics, as we describe below. 

In Refs. [331 EH 1131 [H] (based on the formalism of [15]) results were provided in both 
single-particle-inclusive (1PI) and pair-invariant-mass (PIM) kinematics. The kinematics am- 
biguity was studied in detail in [31] and found to be an important source of uncertainty. In 
1PI kinematics the soft-gluon logarithms are of the form [In fc (s4/m 2 )/s4] + and the soft-gluon 
corrections to the double differential cross section, d 2 a/(dtidu\), are calculated. In PIM kine- 
matics, the soft logarithms are of the form [ln fc (l — z)/(l — z)]+ with z = M 2 / s, and z — > 1 near 
threshold, where M 2 is the it pair mass squared. In PIM kinematics, the soft gluon corrections 
to the double differential cross section, d 2 a /(dM 2 d cos 8), where 9 is the scattering angle in the 
partonic center-of-mass frame, are calculated. The cross section in PIM kinematics was found 
to be smaller than the 1PI result. The results in Refs. [331 El] were based on NLL resummation 
and were later improved by adding subleading terms [131 H3|. The kinematics ambiguity was 
thus reduced in [13^ 114] . Still it was shown in [13j that the PIM kinematics gives large negative 
results at NNLO for the gg channel at LHC energies (for it production at the LHC, the gg 
channel is by far dominant over the qq channel). These negative corrections are deemed un- 
physical since already at NLO the PIM approximation for the corrections does not reproduce 
well the exact NLO result while the 1PI result is a much better approximation (detailed graphs 
for the partonic scaling functions in 1PI and PIM kinematics were shown in [31] and also [T3]). 
In the present paper we have thus used 1PI kinematics. In contrast, Ref. [30] uses a modified 
PIM kinematics. Although the modified PIM kinematics of Ref. [3D] produces less negative 
results than PIM in [31] and [13], the overall NNLO contribution in modified PIM is still neg- 
ative. This explains why both the NNLL resummed cross section and the NNLO approximate 
cross section with modified PIM in [3Q] is less than the NLO cross section at fi = m for both 
Tevatron and LHC energies. This is in sharp contrast to the 1PI results here and in all our pre- 
vious calculations (at both NLL and NNLL accuracy) where the NNLO soft-gluon corrections 
are found to provide a positive enhancement of the NLO cross section. The 1PI kinematics 
provides an excellent approximation as evidenced by the fact that the NLO approximate 1PI 
corrections from the expansion of the resummed cross section account for well over 98% (up 
to 99.8%) of the exact NLO corrections in the gg channel (with fx = m) at both Tevatron and 
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LHC energies. This is a far better agreement than can be attained with PIM or modified PIM 
kinematics. We thus remain of the opinion that the results in [30] do not accurately reflect the 
true contribution of soft-gluon corrections. 

It is also interesting to compare the results in this paper with our previous results in [TSJ [H] . 
Although NNLL resummation requires calculation of the two-loop soft anomalous dimension 
matrices as presented in this paper, it was argued in [T3| IT4"] that the numerical contribution 
of this matrix at two loops to the cross section is expected to be small. In [131 E] many of the 
terms beyond NLL were already included in the calculation and it was argued based on the 
study of the scaling functions in 1PI and PIM kinematics that these additional subleading terms 
were relatively dominant. Now that the full two- loop NNLL terms are known it is important to 
revisit the validity of this argument. We find that indeed the new two- loop terms from the soft 
anomalous dimension matrices contribute very little to the total cross section, and hence the 
argument was valid and the results in [T3"l H4"] were robust. For example, using the MSTW2008 
NNLO pdf [37] the calculation at the accuracy of Ref. [H] for the top quark cross section at the 
LHC at 7 TeV gives 165 pb, while in this paper we find 163 pb based on NNLL resummation. 
The difference between the two numbers is very small compared with the overall theoretical 
uncertainty. Any differences in the numbers provided in [131 H3]> an d the present work are 
overwhelmingly due to the use of different pdf and only to a rather small extent due to the 
different theoretical accuracy. 

To conclude, we have shown in this paper that the top quark cross section and trans- 
verse momentum distribution receive significant enhancements from soft-gluon corrections at 
NNLO. These corrections have been resummed at NNLL accuracy by calculating the two-loop 
soft anomalous dimension matrices for the partonic processes. Approximate NNLO total and 
differential cross sections have been derived from the NNLL resummed result. The NNLO soft- 
gluon corrections enhance the total cross section and the px distribution and greatly reduce the 
theoretical uncertainty from scale variation. The pdf uncertainty of the cross section has also 
been presented. Our NNLL resummation formalism can be used to calculate other differential 
distributions of interest, such as the top quark rapidity distribution. This will be a topic of 
future work. 
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